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Abstract 

The  Watson  transformation  permits  us  to  express  the  radial  dependency 
of  certain  wave  functions  in  terms  of  Hankel  functions  of  complex  order. 
These  Hankel  functions  may  be  considered,  as  "eigenfunctions"  of  a  non-self- 
adjoint  boundary  value  problem  for  Bessel's  differential  operator.   It  turns 
out  that  in  general  an  arbitrary  function  cannot  be  expemded  in  a  series  in- 
volving these  eigenfunctions,  not  even  if  very  strong  conditions  are  imposed 
on  the  arbitrary  function.  Hovever,  there  exist  sufficiently  'weak'  non- 
selfadjoint  perturbations  of  a  selfadjoint  problem  connected  with  Bessel's 
operator  for  which  an  expansion  theorem  for  an  arbitrary  function  in  tenns 
of  the  eigenfunctions  can  be  proven. 
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1.   Introduction 

In  the  theory  of  propagation  of  electromagnetic  waves  aroiind  the 
surface  of  a  sphere  produced  hy  a  unit  dlpole,  the  field  can  he  expressed 
in  a  series  of  Hajikel  fiancfions  H    (kr).   These  are  of  Integral  order  n) 
the  variable  is  kr  where  k  is  the  wave  number  and  r  is  the  distance  from 
the  center  of  the  sphere.   If  kr  is  large^  this  series  converges  very 
slowly  and  is  tinsultahle  for  numerical  computations.   It  has  been  shown 
by  G.N.  Watson  and  others  that  in  this  case  a  convergent  expansion  for  the 
wave  (at  least  in  the  "shadow  region")  can  be  obtained  in  the  form  of  an 
infinite  series: 


oo 


u  =  y]  c  H  ^^^(kr) 
^^   n  V   ^ 


n=l 


where  the  numbers  v  are  the  roots  of  a  transcendental  function  of  v.   In 

n 

the  case  of  a  perfectly  conducting  sphere,  this  function  is  simply 

(1.1)  H^^^)(ka). 

It  is  known  that  there  exists  an  infinitude  of  complex  roots  of 
H  (ka)  =  0. 

V 

It  has  been  pointed  out  by  Sommerfeld  ^  -I  that  the  functions  H  (kr) 

n 
form  an  orthogonal  system.   If  we  replace  kr  by  a  real  variable  x,  and 

assume  that  k  =  1,  we  have 

oo 

(1.2)  f    H^(l)(x)H;^)(x)f  =0 

^     n       m 
a 


if  m  5^  n. 


-   2 


Therefore,    if  a  function  f(x)    can  be  written  in  the  form 

oo 
(1.5)  f(x)   =     ]C     ^n^^^^^^) 


n=l  n 


and  if  the   series  on  the   right-hand  side   converges  absolutely  and  uniformly, 

ve   can  determine  the   constants   c     from 

n 

00 


^  n 


{l.k) 


a 


n  oo 


/   «^')  - 


[21 
B.  Friedman  '-  -■  studied  the  whole  complex  of  questions  arising  from 

the  Watson  transformation  '-  -^ .  He  showed  that  the  question  of  whether  f(x) 
can  be  expanded  in  a  series  of  the  type  (1.5)  inay  he  considered  as  a  gener- 
alized (non-self ad Joint)  eigenvalue  problem  for  the  second  order  differen- 
tial operator  L  defined  by 

/n  =^  T     2  d^     d     2 

(1.5)         L=x  —2  +  ^0^  +  ^ 

dx 

where  the  boundary  conditions  are  given  for  the  interval   (a, oo).      The   con- 
dition at  X  =   a  is  aji  ordinary  boimdary  condition,   but  the   condition  at 
Infinity  is  of  the  type  of  a  "radiation  condition"    (see   (I.7)).      B.   Fried- 
man suggested  that  the  problem  of  expanding  f(x)   in  a  series   of  "eigenfunc- 

tions"     H^         (x)   of  this  bovmdary  value  problem  be   treated  by  ujsing  the 

'''  [51 

methods  employed  by  Tltchmarsh   '-  ■' .     This  method  hinges  on  the  proof  that  a 

certain  integral  in  the   complex  v-plane,    taken  over  a  large  semi-circle. 
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tends  toward  zero  as  the  radius  of  the  circle  tends  toward  infinltyo   It 
is  somewhat  siirprising  that  this  method  does  not  work  in  the  present  case. 
This  fact  can  be  demonstrated  indirectly  by  showing  that  the  integral  in 
the  complex  v-plane  will  not  tend  toward  zero  and,  directly,  by  showing 
that  the  series  (I.5)  with  coefficients  (lA)  will  not  converge  for  a 
"reasonable"  function  f(x)  which  satisfies  the  proper  boundary  conditions 
at  X  =  a  and  x  =  00  .   So  far,  no  set  of  necessary  and  sufficient  condi- 
tions is  known  which  could  be  iised  for  deciding  the  question  of  whether  a 
function  f(x)  admits  an  expeinsion  (1.5)  (excepting,  of  course,  the  tri- 
vial condition  that  f(x)  is  given  as  a  series  of  the  type  (l.j)). 

Under  these  circumstances  it  seems  to  be  of  some  interest  to  find  out 
whether  there  exist  any  non-selfadjoint  boundary  value  problems  connected 
with  the  operator  L  (and  of  the  type  studied  by  B.  Friedman)  to  which 
Tltchmarsh's  method  can  be  applied.   This  will  be  done  in  Section  2,  where 
a  dependency  of  the  operator  L  on  the  eigenvalue  parameter  is  assiamed  in 
such  a  manner  that  the  problem  may  be  considered  as  a  perturbed  Sturm- 
Li  ouville  Problem.   (An  ordinary  St\irm-Ll ouville  Problem  connected  with  L 
arises  If  we  consider  the  Interval  (la,  loo)  on  the  imaginary  axis;  in 
this  case  the  eigenvalues  lie  on  the  Imaginary  v-axls).   In  the  case  of 
the  pertiirbed  problem  the  eigenvalues  v  *  may  have  arbitrarily  large  real 
parts,  but  they  will  not  move  as  far  away  from  the  imaginary  axis  as  the 

zeros  V  of  the  function  R  ^   (a)  for  a  real  value  of  a. 
n  V   ^  — 

In  Appendix  II,  we  shall  show  that  the  series  (1.5)  with  coefficients 
(1.^)  diverges  for 


hJ^^(x)   H^^^^(x) 

(1.6)       f^(x)  =  -P -TIT— 

H  ^^^(a)   H  ^^^(a) 


vhere   p  and  (Tare  different  positive  numbers.      This   is   true  although  the 

function  f  (x)   vanishes   for  x  =  &,    is  analytic  and  satisfies  the  radiation 
o 

condition 


(1.7) 


11m 
X  — >  oo 


df  (x) 

—^ if  (x) 

dx  o^ 


X  =  0. 


2.        A    'Perturbed'    Stujnn-Llouville  Problem 


Theorem  1. 


If  f(x)  has  a  continuoiis  second  derivative  in  (a^N)  such  that  f(x)  =  0 
for  x  =  a,  a  >  0  and  for  x  g  N  then  in  the  interval  (a,N)   f(x)  may  be  writ- 


ten 


f(x)  =  Y2     c(v/)?^Jx), 


d   (x)  are  the  eigenfunctions  and  the  v   are  the  eigenvalues  of  the  dlffer- 


ential  operator: 


i.e. 


L  =  x  — g  +  X  —  -  X  [l  +  g(v)] 
dx 


U   (x)  =  V  U   (x) 
^n^      n  ^n\ 


vlth  the  boimdary  conditions 


^J^)   =  0,        ^J^x) 


la  X 

n 


V^ 


as  X  — >  oo  . 


a  =  a(v)   is  defined  by  a  =  ±/T~+~g(vJ,     Im  a  ^  0.     g(v)  is  assumed  to  be 
an  analytic  function  of  v  in  Re  v  >  0  such  that  Re   g(v)   +  1  >  b   >  0 


I  g(  V )  1    g  M       as    I  V  I   — >  00 , 


Im  g(v)  =  er 


dv  ^  V 


log  r  log(log  r) 
r 


and  <  0 


> 


If 
V  =   re        ,    €   >  0 


;(iT)   =  g(-lT)    for  real  t. 


See  Appendix  I, 


The   eigenfunctions  except  for  a  multiplicative  function  of  v  will  be 

H  ^^\a  x) 
V       ^   n   ' 
n 

i.e.  Hankel  functions  of  the  first  kind  vith  v  chosen  so  that 

n 

n 
We  will  omit  the  superscript  (l)  in  vhat  follovs  and  vrite  simply  H  (ax). 
Proof: 

Consider  two  auxiliary  functions  "^^{x)   and  "i^^x)    such  that 

(L  -  v^)Y^(x)  =  (L  -  v^n^ix)   =  d 

T  (a)  =  0       (to  satisfy  the  lover  boimdary  condition) 
Y_(x)  =  H  (ox)   (to  satisfy  the  upper  boundary  condition)  . 
We  choose  f, (x)  as 

Y^(x)  =  J_^((xa)j^(oQ£)  -  J^(aa)J_^(ax). 


Consider  the  inhomogeneous  eq-uatlon 


(L  -   v^)  f{x,v)  =  f(x) 


where  f(x)  satisfies  the  conditions  of  Theorem  1. 


-  6 


J(x,v)  may  "be  expressed  in  terms   of  "f^Cx)   and  ^2(^)- 

-jtH  (ax)  p      J     (aa)j   (ay)    -  J   (aa)j     (ay) 


a 


rtJ_^(aa)j^(aoc)    -  J^(aa) J_^(cuc)        p        H^(ay)f(y) 


2H   (aa)slnjtv  J  y 


/ 


Consider 


^       /       $(x,v)v  dv 


C 
m 

where   C     is  an  appropriately  chosen  sequence   of  contours  in  the  right  half 
plane.      (We  may  consider  them  as  a  sequence  of  approximate   semi-circles, 
the  radii   of  vhich  tend     toward  infinity. ) 

We  observe  that  along  the  imaginary  v-axis      v^(x,  v)    is  an  odd  func- 
tion of  V)    hence 

iT 
~  \     /  $(x,  v)vdv  =  0.  (See  App.   1,    Sec.   Ill  and  IV^ 

-iT 

^(x, v)  is  a  meromorphic  function  of  v  in  the  right  half -plane.  .Hence 
we  may  evaluate  the  Integrals  hy  considering  the  residues  due  to  the  zeros 
V  of  H  (aa).  At  these  points  the  wronskian  of  '^-.(x)  and  Ifp(x)  vanishes. 
Therefore 

f^(x)  =  k  ^^(x) 

where  k  does  not  depend  on  x. 
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Hence 


^  /  $(x,v)vdv  = 


/  I  c^(v^)H^^  (ox)    (a  finite  sum), 


n  n  V 


n 


m 


On  the  other  hand,  we  will  now  show  by  considering  the  asymptotic  form 
of  ([(x,  v)  for  large  |v|  that  the  above  integral  over  the  sequence  of  contours 
C  tends  toward  f(x)  as  m  — >  oo  .   For  this  purpose  we  will  need  the  follow- 
ing: 
Lemma: 

For  large  |v|  the  zeros  of  H  (aa)  in  the  right  half  plane  have  their 
argument  asymptotic  to 


2 


h(v) 


1  +  lOf 


2r 
Aea 


where 


a  =  Ae 


19 


V  =  i^ 


If 


Proof: 


9=1  -  h(v). 


For  I  V I  — >  oo 


H  (aa)  = 


J  (<xa)  -  J  (aa)e 

—  V  V 

1  sin  nv 


-jtvi 


1  sin  jtv 


(f) 


/oax    -Jtvi 
(-)   e 


_r(i-v)     r(i+v) 


and  therefore  if  H  (aa)  =  0, 

V 


-  8 


-2v  „    „      ^-rtvi 
,aB.^  2v     -2v     e 

(-5-    V   e    -^ 


2  sin  jtv 


Taking  absolute  values  we  see  that 


,Aea  "^''''°^*  2(9-Y)rsinV? 

I )        e  ' 

^2r  ^ 


e    -  1 


Since  1.'''°°^  dominates  for  r  — >  00  in  most  of  the  upper  right  half 

plane^  it  is  clear  that  cos  1   must  be  small,  or  that  "i     =  arg  v  when 
H  (aa)  =  0  is  of  the  fonn 

^0=1-^ 

where  5  is  small.  Hence  we  may  replace  cos  "F  by  6  and  sin  T  by  1  and 


we  find 


or 


2r 
6  log  1^  ♦ 
°  Aea 


^o~2 


[|-Mv)-|.5]->0 


h(v) 


1  +  log  ( 


Aea^ 


(cf.  App.  1,    Sec.  IV.) 


Note  that  to  the  right  of  the  ray  arg  v  =  1  ,    J_  (aa)  is  the  dominant 
part  of  i  sin  nv   H  (aa),  and  between  the  positive  imaginary  axis  and  the 
ray  Y  =  f  ,  J  (aa)e~    is  the  dominant  part. 

This  leads  i:is  to  consider  a  similar  ray  in  the  v-plane  on  which 


J  (2Na) 


J  (2Na)  e 


-rtvi 


We  have  for  ] v [  — >  00 


9  - 


J  (ccx)e 


■rtvi 


J     (ax) 


2v     „, 

(f)     ^ 


,CQCv 


2v 


-rtvl 


vr^( 


vjsinjtv 


2rcos^     2r(f-9)8in* 

^2r   -^  2itvl      , 

e  -   1 


X 

2W 


2rcosf      „    ^,,  2rco3T     2r(^-9)sinY 
^    2r   '^  I    2jtvl      ,  I 


^   (— ) 


2rco8'] 


M 


(cf.  App.   I,    Sec.   I.) 


vhere  M  is  independent  of  both  v  and  x.  '  • 

This  is  valid  in  the  region  between  the  ray  described  above  and  the 
imaginary  axis,  assuming  that  v  remains  a  bounded  distance  avay  from  the 
integers.  Hence  for  any  e  >  0 


_  /   V  -rtvi 
J  (ax)e 

V 

J  (ax) 

-V 


for  all  X 


<  e 


a  S  X  S  N 


on  a  contour  r  vhich  ve  vill  take  as  a  seml-clrcle  of  non-integral  radius 
extending  from  the  negative  imaginary  axis  to  the  line  vhere 


J_^(2Na) 


e-^^^J  (2Ka) 


(sufficiently  large  v) 


10  - 


This  simplifies  our  asymptotic  expression  for  $(x, v)  since  ve  may- 
replace  the  Hankel  functions  by  tlie  simpler  Bessel  functions  J^  on  P^. 
TUmB  ve  see  tliat 


l(x,v) 


-rtJ      (ox) 
-v" 

2J      (ixa)sinjtv 


jt  ('j_^(cLa)j^(cuc)    -  J^(aa)j_^(a^)j 


2J     (aa)    sin  rtv 


N 


J_^(ay)f(y)dy 


A. 


Y\V 


^x 


(-) 
^xy' 


f(y) 


dy 


N 


-/ 


(-)      -   (— ) 
y  ^xy' 


f(y) 


dy 


Consider 


/ 


(Z)     £(XI  dy  =  i(I)    f(y) 


y=x 


y=a 


Integrating  the  last  term  by  parts,   ve   see- that 


/    (J)'^dy  =  if(x).-^(|vr2) 


We  note  that  the  other  terms  of  the  asymptotic  expansions  of  the 
Bessel  functions  do  not  contribute  a  lower  order  tenn: 


-  n 


21  V 

(^)  (^)   /      2  2 


Applying  a  similar  argument  of  integrating  by  parts  to  th.e  terms 
resulting  from  consideration  of  these  additional  terms^  we  conclude  that 
no  —  term,  i.e.  larger  term,  enters.   Similarly, 


/  ¥  ^  ■  I .... .  it^) 


and 


N  N 


Hence 


-  il  jrr^$(x,v)vdv~-^i^  f(x)  — >f(x) 

as  the  radliis  of  the  sequence  of  semi-circles  tends  towards  infinity.   6' 
is  the  angle  that  the  ray  ■    .  •.  u  i, 

I  J   (2Wa)  I   =   I  J  (2Wa)e~^^^ 
makes  with  the  imaginary  axis. 


We  now  choose  contour  Pp  so  that 


arg  J  (aa)  -  arg  [  J  (aa)e     j  =  const  4   2nl 


where  I  is  any  integer. 


12 


This  will  insure  the  denominator  being  bounded  from  below  as  the  con- 
tour Pp  passes  through  the  region  near  the  zeros  of  H  (aa).  Investigating 
this  further,  we  find  that  If 


V  =  cr  +  It 


then 


2v 


^^g  e-«^i  (^)   ^  ^  2(59  +  a-r  log  ^  -  T  log(cy^  t^) 


¥J^ 


2cr  tan"  ~+   2t  +  ^. 


(We  observed  that 

arg  r(v)  =  Im  log  r(v) 
and  applied  Stirling's  form-ula. 


*"■      -  log(<r"+  .")  +  log{f )  .  20-g  +  m 


This  is  small .   Hence  r^   is  approximately  a  straight  line,  the  length  of 
which  is 

^  const  r6'  ~  const  rh(v)  ^ 

Along  that  part  of  P^  to  the  right  of  the  ray  where 


15 


ve   can  find  a  constant   so  that 


J      (aa)    -  J   (aa)e 


-nvi 


^   const 


J     (aa) 


This  ve  will  designate   as   P- 


The 


remaining  part  of  V^,    denoted  by  T,  ,   will  have  the  property: 


J      (aa)    -  J   (aa)e 


-nvi 


a    const  J   (aa)e 


■jtvi 


Considering  as  before 


J   (ax)e 


-Jtvi 


~J     (ax) 


we  conclude  on  Tr : 


H  (ax) 
v___ 

H  (aa) 
v^   ' 


^  const 


J  (ax) 

J  (aa) 
v^   ' 


rS 


^  const(— ) 
^a 


On  r. 


H  (ox) 

V 

H  (aa) 

V 


S  max   const 


J   (ox) 
-v^ 

J   (aa) 

-V 


,  const 


J  (ax)e 

V 

~J  (oa) 


-Jtvi 


^  majc 


-Rev  Rev 

const  (— )     ,    const  (— ) 


r6' 


g  const  (— ) 
^a 


Ik 


since 


J  (ax)e 

V 

J   (aa) 

-V 


g  const 


J  (ax)e 


-itvi 


J  (aa)e 


-Ttvi 


on  r. 


Consider 


(f) 


(-) 

y 


^  dy  =  i  (f )   f(x)  +  i  (|)   f(x)  +  cr 


M^^^ 


V  a 


V   X 


=  ^|  -i-Ml^l^' 


vhere  M  Is  some  constant.  Hence  if  w(y,  v)  denotes  the  vronskian  of  *-,(y,  v) 


and  T„(y, v),  then 


J  rgVdvfgCx.v)  J 


*l(y.v)f(y) 
w(y,v)y 


<ay 


S   const    I  V I  .  I  V 1    6 '  — ^  M^^ ' 


v6' 

=  const  S'M 


This  vill  tend  to  zero  as    |v|    tends  toward  infinity  if 

log(log    |v|) 


6'  =  cr 


since 


l^lB,  /  log(log    |vl)  \ 

=  ^(   -j j   (log    |vl)l°SM^ 


const  6'm'   '      =  er 


15 


This  vill  be  true  vhen  g(v)  satisfies  the  conditions  of  the  theorem: 


a  (v)  =  -  (l  +  g(v) 


2  210 
A  e 


For  large  | v ]  ve  have 


^   1    4-     Im  g(v)     n 

9  ^  y;   arctan  ;; /■  V  ■  ■  ■  +  tt. 

2        Re  g(v)  +  1   2 


(We  assume  that  IM  g(v)  is  negative  in  this  region.)   Since  Im  g(v)  is 
small  and  Re  g(v)  is  bounded  avay  from  -1,  we  have 


-h(v)  +2=^~2Reg(vj  -;i^2- 


(.: 


Therefore, 


h(v)     =  0-  I  log  ^  log(log  r)  ^_  ^ 
1   +  log  T^     \  r  f  1  +  lo. 


AeW 


(l  +  ^°^A^) 


log (log  r) 


We  observe  that  this  Inrplies  that  the  real  parts  of  the  v  may  tend  toward 
infinity. 

A  similar  consideration  is  valid  for 


N 


n  n      ^p(y,v)f(y) 

/     vdv.Y^(x,v)   / 


Ky;.v)y 


dy. 


Therefore  we  conclude  that 


11m    ^   /   5(x,v)vdv  =  f(x) 


m  — >  oo 


m 


-  16  - 

Evaluating  the  Integral  "by  the  theorem  of  residues  we  have  therefore 

f(x)  =  J2     C(v^)  H^  (cue) 

n 

for  any  x  in  the  interval  (a,N). 

(The  zeros  of  sin  jtv  do  not  make  any  contribution  to  the  integral 
since  at  the  integers 

J  (z)  =  (-1)''  J  (z). 
n^     ^      -n^ 

Hence  at  these  points   the  niamerator  simultaneously  Vanishes. ) 

Theorem  2. 

If 

~  2 

L  =  L  -  q(x)x 

vhere  L  is  the  operator  of  Theorem  1  and  q.(x)  is  a  continuous  function  of  x 
in  (a^T)  such  that 

q(x)  =0  X  g  T; 

then  any  function  f(x)  satisfying  conditions  of  Theorem  1  may  be  expanded 
in  terms  of  the  eigenfunctions  of  the  operator  L  with  the  same  boundary 
conditions. 
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Proof; 

We  will  denote  the  flmctions  of  Theorem  2  that  correspond  to  those  of 
Theorem  1  by  a  tilde « 

Consider  *^(x, v)  such  that 

(L  -  v^)l^ix,v)   =  0 

^^{a,v)    =   0. 

^, (x, v)  may  be  written  in  terms  of  ^   {x,v),    the  auxiliary  function  used 
in  the  proof  of  Theorem  1: 

X 

f^(x,v)  =  f^(x,v)  +2li^  J    y(j^(a-x)J_^(a-y)  -  J_^,(c^)J^(ccy))g(y)i^(y,v)dy, 

a 

For  large  |  v  | 


V       V 


X    I- 


V^^^)~  r(i+v)r(i-v)  • 


1  +  & 


*■  i  I   y 


V       V 


i.^a) 


q(y)^i(y.v)dy 


-V 


(f)  v^'^) 


lr(i+v)r(i-v) 


X 

J    y  (J)  fj^(y,v) 


k(y)  l^iy. 
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Let 


M  =    max 

a  g  y  ^  X 


Then 


\ 


(J)  ^^(y.v) 


^M^Q 


|r(i+v)r(i-v)|   Ivl 


for  sufficiently  large  \v\,   vhere 


Q  =  J  yl(i(y)l<iy. 


Hence 


|r(i+v)r(i-v) 


"i-  — ^ 


8 


1  - 


|r(n-v)r(i-v) 


assuming  | v |  S  8q. 

SutBtituting  this  bound,  we  find  that 


const 


(f) 


?^(x,v)  -  Y^(x,v) 


|r(i+v)r(i-v)| 


Hence,  except  in  the  vicinity  of  the  imaginary  axis,  ve  have 


?^(x,v)  =  f^(x,v)  (  1  +  e^( 


Similarly,   we   consider 


f2(x,v)   =  ^2^^'^) 


2sinrt 


X 

^    J     [j^(a,x)j_^(ay)    -  J^(ay)  J_^(ax)1yq(y)l2(y.  v)dy. 
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We   consider  first   v  along  P      as   defined  in  the  proof  of  Theorem  1,   vith  N  re- 
placed by  max   (N,T). 


?2(^^^) 


{^)     (l-K>(l))- 
ir(l-v)    sin  nv 


T   r 


1 

V 


X 


V  V 

(f)  -  (?) 


'X 


1  -^  CT(i) 


yq(y)^2^^'^^'^"^' 


(^)  ?2(x,v) 


r(l-v)sin3tv 


(^)  ?'2(y,v) 


•y((i(y)|dy. 


Let 


M^  = 


max 
X  ^  y  ^  T 


(f )  yy.v) 


Then 


M                          2                        ^^2 
M     S  +  

|r(i-v)    sin  nv|         jvl 


Hence 


M     ^   ir'(l-v)    sin  rtv 
2  -  ^  _   4Q_ 


|r(l-v)   sin  rtv| 


for    |v|    ^  8Q  as  before. 

Substituting  this  bound,   ve. obtain: 


fgCx^v)    -  ^^{x,v)\    ^ 


const 


l(¥) 


|v|  |r(l-v)   sin  jtv| 
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Hence  along  r^  except  in  the  vicinity  of  the  imaginary  axis^  we  have 


l^ix,v)   =  'J'gCx^v)  (  1  + 


Along  r  ve  have: 


(^)\(x,v) 


(^) 


sin  jtv 


max 
X  g  y  S  T 


(f) 


-V 


r(i-v) 


(f )  e-'^^i 


r(i+v) 


-  f 


(f )  ^^ij,y) 


yk(y)  Ky" 


Hence  along  P^ 


^2(^^^) 


g 

4 

(f)^ 

sJ 

-n  nv 

max 
X  S  y  S  T 


(^) 


-V 


r(i-v) 


(f )  e-'^^i 


r(i+v) 


We  nov  obtain  estimates  for  the  wronskian,  of  (^f  (x^v)¥p(x,  v)j  which 
will  be  denoted  by  c3(x,v). 


d?  (x,v)   ^       d?  (x,v) 
;:(x,v)  =  Y^(x,v)    ^    -  V^'^)    dx 


=  W  (j^(x,v),T2(x,v)) 


+  W 


1',  (X,V)  -  ■;r ^ 

1^  '     2  sin  Jtv 


G(x,y)>i'2(y,v)yq(y)dy 


-    21 


+  W 


g  3^"^  ^^     /     G(x,y)?^(y,v)yq(y)dy,    1^{y,v) 


+  W 


2  sin 


.A. 

r^^^    J      G(x,y)i^^(y,v)yq(y)dy, 


2  s 


^^^   J  G(x,y)?2(y.v)yq(y)dy 


=  W^  +  W^  +  W     +  Wj^, 


in  which  w[]      j    denotes  the  wronskian  and 


G(x,y)   =  J^(cQc)j_^(ay)    -  J^(ay) J_^(ax) , 


We   observe  that 

W^  =  \i{l^{^,v),lj,^,v))   =  cj(x,v) 

is  the  Wronskian  of  Theorem  1. 
We  vlll  now  show  that 


w^ 

W 

w, 

2 

> 

and 

1 

1 

1 

are  small  along  both  P.,  and  Tp  and  therefore  we  will  be  later  justified  in 
replacing  w  by  w. 

We  observe  that 
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Wg  =  w 


T 


hi^'^")'  2  sTnnv   f  G(x,y)?2(y^^)y<i(y)<iy 


2  sin 


T 


+  -  /   G(a,y)"iJy,v)yq.(y)dy. 


Along  r,  for  sufficiently  large  |v|: 


const 


/ 


(i) 


(^)  yk(y)|dy 


r(i+v)r(i-v) I   |r(i-v)  sin  rtv| 


/OcLs 

\  2  > 


Ir(i-v) 


const 


Along  r,  we  have: 


T     Ly 


const 


(t) 


!(|)  ||(f)  I  yk(y)|dy 


lr(i+v)r(i-v)|     lr(i+v)  sin  ^vj 


r(i+v) 


const 
g  M     where  M  Is  some  constant. 


M^"^ 


|v| 
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We  conclude  on  T,  ^  since  Re  v  =  erflog(log  r)j  ,  that 


->  0     as   V  >  00 


Along  r,  ve  have 


const 


(^)    (-) 
^a    V 


max 


W, 


W, 


|r(l+v)r(l-v) 1   y  s  z  g  T 


(f) 


-V 


r(i-v) 


(-5-)  e 


r(i+v) 


yk(y)|<iy 

I  sin  -nv  \ 


/Oax 
^  2  ^ 


r(i-v) 


const 
g  M     as  in  the  previous  case. 


Therefore J  we  conclude  that 


W. 


W, 


>  0  along  the  entire  contour  C  as  m  — >  00 ■ 

m 


We  treat  the  ratio 


W. 


W, 


in  a  similar  manner: 


W  =  W 


2   sin  Kv  /  G(x,y)?^(y,v)yq(y)dy,'i'2(x,v) 


.A. 

2  sin  ^v  /  w[G(^^y)^^2^'''''^]^l'^y''')y^^y)'^y 


2h  - 


1 

X 


H^(ay)?^(y,v)y(i(y)dy. 


We  have  along  T  for  sufficiently  large  |v| 


const 


W, 


W, 


^      lr(i 


(^) 


(f)Vk(y)|dy 


|r(i-v)siri  nv|   |r(i+v)r(i-v) 


I   -V 
(f) 


r(i-v) 


const 


Along  r  and  V^^   we  find  as  in  the  previous  case 


W. 


W, 


const 
S  M       for  some  constant  M. 


W, 


^       k  sin  nv 


rt 


W 


X  T 

r    G(x,y)?j_(y,v)yq(y)dy,    r    G(x,y)?2(y.  v)yci(y)dy 


2x  sin  jtv 


/ 


J^(ay)f2(y,v)yq(y)dy 


J      J_^(ay)f^(y,v)yq(y)dy 


f      /     J^(ay)?j_(y,v)q(y)ydy  J    .    f     /     J_^(ay)f2(y.  v)q(y)ydy  j 
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Substituting  the  bounds  as  done  previously  we  find: 


\ 


^1 


const 

g  in  r, 

lv|2  1 


const  _, 
^  P  M^"'  in  r. 


->  0  along  C  as  m  — >  oo  , 

m 


Therefore,  we  conclude  U)(x,  v)  =  W  +  \1  +   W  +  W,  is  dominated  by  W  =  w(x,  v) 


and  along  the  entire  eoutour  C  and  specifically: 


to  =  u)(  1  +  er(  1  I        along  P 


'M^^^ 


CO  =  CO  I  1  +  Cf  I  1  1        along  Tg. 


Consider  the  inhomogeneou^  equation 


(L  -  v^)f{x,v)   =  f(x). 
^(x, v)  may  be  expressed  in  terms  of  ^^(x, v)  and  Yp(x, v): 

r  ^-,iy>^)^iy)^y    ^  r    ^2(y>v)f(y)dy 

J(x,v)  =T2^x,v)   /  -— _^_  +  f^(x,v)   /   — — -^- 

^^    ^y;.v)y  ^      u)(y,v)y 

We  will  consider  this  expression  for  large  v  and  show  that: 


11m    ^  /  |(x,v)vdv  =  f(x), 
— >  00  ^'^    ^^ 


m   .  ^     C 

m 
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Consider 


"S^ix^v) 


¥^(y,v)f(y)dy 
w(y,v)y 


''^2^'''^^    "   2  sin  jtv    /    G(x,y)?2(y.v)yq.(y)dy 


^1^^'^)   +  2  sin  nv    /     G(y.z)i^(z,v)zg(z)dz 


f(y)'iy 

J(y;>v)y'' 


=  Y  (x,v)  r  4 

^  (.1 


^l(y,v)f(y) 


X 


'^(yjv)y' 


2  2  sin  rtv     2 


^o(^^^) 


f(y)dy 

r 

'^(y^v)y'' 


X      /      G(y,z)i^(z,v)zq(z)dz 


2  sin  jtv 


G(x,y)^„(y,v)yq(y)dy 


'i-^(y,v)fiy)dy 


X 


wl 


(y;.v)y'' 


+  - 


h  sin^jtv 


f    G(x,y)i   (y,v)yq(y)dy  .  r      ^^^^^      T    G(y,  z)?  (  z,  v)zq(z)dz 


a  '^(y'^)y    a 


=  |^+$2+$5+li,    • 
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We  note  that 


r     ^.(y,v)f(y) 
J^=^^(x,v)      /   -i- -5-  dy 


"^(y.  v)y 


.A. 


^l(y,v)f(y) 
u(y^v)y 


dy 


which  was  considered  in  proof  of  Theorem  1.   Hence 


lim    — r  /  ^  ^  vdv  = 


m  — >  00 


f(x) 
2 


m 


We  will  now  show  that  the  contributions  to  the  contour  integral  due  to 
Jp,  ^   and  ^  ,  tend  toward  zero.   Consider 


X  y 

I2  =F-iI^  V-'^)  /   rr^  /  G(y,.z)?^(z,v)z<i(z)dz. 


a    '^(y'^)y    a 


We  interchange  the  order  of  integration  and  note  that 


yw(y. v) 


is  a  function  of  v  only. 


I 


2        2  sin  jtv 


yw(y,v) 


?^(z,v)zci(z)dz  I    ^iy'f(y)^y 


Consider: 


.A. 

/ 


X    ,- 


G(y,z)-f(y)dy 

y 


r(i+v)r(i-v) 


V  V 

(?)  -  (p 


f(y)dy 

y 


sin  Jtv 


Jtv 


V  V 

(^)    +   (-) 


,X  X     _ 


f(y) 


/ 


V  V 

(^)  +  (-) 

^z'        ^y' 


f'(y)dy 
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by  integrating  by  parts.      Hence 


G(y,z)f(y)dy 

y 


const     sin  rtv 


(f) 


Therefore  along  P^ 


const 


(?) 


(i' 


% 


Isin  jtvl         |r(l-v)sin  nv|       ^        |r(l+v)r(l-v; 

9, 


const 


(f) 


|r(i-v) 


sm  nv 


(f) 


dz 


5   ' 


Along  r,  we  have 
5 


ft 


const 

max 

(^)" 

V 

(2)    - 

V 

sin  jtv 

sin  rtv 

r(i-v) 

r(i+v) 

•          v' 

-V 

,aa\ 
.2  ^ 

|r(i-v) 


M^" 


g   const  =■     for  some   constant  M. 


In  like  manner  we  obtain  this  bound  for    ^1    along  Tl  . 


-  29 


Hence 


^  vdv 


const 


/  f^^dv 


S  const 


M^ 


ev 


and  we  conclude  that 


11m     /  J  vdv  =  0. 


m  — >  00 


m 


Similarly  ve  find 


fel 


2  2 


h   sin  jtv 


G(x,y)5'Jy,v)yq(y)dy 


r   f(y)dy   " 


^a^y.v)y   ^ 


2  /  G(y,z)Y^(z,v)zq(z)dz 


const 


for  V  along  r^ 


const  ,J{ev 


M^^ 


for  V  along  r^ 


and 


15; 


2  sin  3tv 


G(x,y)T,(y,v)yq(y)dy 


/ 


\(y,v)f(y)dy 
'^(y^-vjy 


const 


for  V  along  P, 


const  .J^ev 


5 


M«^ 


for  V  along  r 


2* 
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We  conclude  therefore 


lim    /   (  ^,  +  ^)vdv  =  0. 
-  -^«^  Cm 


In  the   same  manner  we  find 


11m 


N 

TtX  1  ,/  _.  n  -  o 


m  — >  CO  -^^  /^2  2 

Cm  X  '^{y,v)y 


and  therefore 


lim     ^-T-   /  ^  (x^v)vdv  =  f(x). 
m  >   00      "-^ 


As  in  Theorem  1^  we  also  evaluate  this  integral  by  the  theorem  of 
residues.   Hence  we  conclude  for  x  in  (a,N)  that: 


f(x)  =^  c(v^)  ?^(x) 
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Appendix  I 


In  this  appendix  modifications  of  the  hypotheses  of  Theorem  1  will  be 
discussed. 
I.   Let  us  consider  the  following  conditions: 


1) 


a(v)  =  ef(v^-^) 


for  some  positive  h  as  |v[  — >•  oo 


replacing  hoimdedness. 

2)   Instead  of  a  wedge  we  consider  the  region 


2  -     M    -  2 

k 

1-c 
log  r 

and 

3) 

a'(v)      =  o 

1,,                     vl-C 

-(log  r) 

for  some  positive  constants 
k  and  c,    as  |v|  — >  cD 
in  this  region  for  the  same  constant  c. 

Then  we  note  that  the  asymptotic  formulas  for  the  Bessel  functions  are 
still  valid  for  example: 


J  (ax)  = 


(f) 

r(i+v) 

2  2 
-,          ax, 

•   • 

(f)' 

1  +  Ot    V  -^) 

r(i+v) 

• 

Furthermore  it  is  necessary  to  show  over  F^  that 


52 


<  €. 


2rcosT  „  ^,  2rcos'y   2r(?-9)sin^ 
_  (iL)       (A£N)      .  e_; ^ 


This  is  clearly  small  in  most  of  the  right  half  plaae  but  near  5 ' , 
it  may  not  be.   In  the  vicinity  of  6'  this  expression  is 


2rcosf 


(— ) 


exp 


2r(-  6  log  -^  -  5  +  h(v)) 


Consider 

f(6)  =  -  6  log;^  -  5  +  h(v). 

We  require  that  f(6)  <  0  for  6  >  5'   [f(6')  =  0,  6'  same  as  in  Theorem  1.] 
If 


5  S 


^l-C 


(log  r)' 

then  clearly  f(6)  <  0  for  sufficiently  large  r.  For  6  ^ 
We  compute  the  derivative 


d6-  =  -  ^°S  J^jj  +  6  -  ^  -  1  +  ^p. 


(log  r) 


1-c  • 


and  vlll  show  that 


55  - 


df(5) 


d6 


<  0 


in  this   range. 
Consider 
a(v) 


A( v)exp 


•  /rt 


i{%-   h(v)) 


Sa 
^ 


iv  = 


1  Sa 

A  36 


,  Mid' 


Hence 


1  bA 
A  'SE 


=   -  Re 


[^] 


ah    ^ 


Iva'l 
.    ^  J 


We    conclude   that 


d5 


f(6)   <  0 


when 


|a'(v) 


)  -(log  r) 


We  note  that  this  bound  for  |a'|  also  suffices  to  make 


d(r 


small  so  that  we  may  pass  a  curve   P    between  the   zeros   and  obtain  proper 
estimates. 


II.     We  may  relajc  the   condition  that  Im  a(v)   ^  0  in  the  right  half  plane 
by  merely  reqTiiring  that 


Im  a(v   )   SO. 

^   n 


-  3^ 


Furthermore  should  a(v)  have  an  infinite  singularity  on  the  Imaginary  . 
axis  then  ve  require  the  exigence  of  a  sequence  of  curves  converging  toward 
the  point  and  on  which  Im  a  ^  0.      To  Justify  this  we  note 


yx)   =  H^(ax)  ^/^  exp[l(ax  -  f  -  5)_ 

f    (x)   =  J     (aa)j   (ox)    -  J  (aa)j      (ax)  ~  -|     / -i-  sin  jtv  sin  a(x-a) 


for  large  a. 


Consider 


^    f,(y,v)f(y)                                 f  ^p(y,v)f(y) 
$(x,v)   =,^J'2(x,v)      /       — 5-dy  +  f^(x,v)      /  ^—  dy 

'-'  i.il'-u-.vW  '^  (,/v.vW 


e 


w(y,v)y  '-'^       w(y,v) 

la(x-a) 


V 

x 
X  N 


r    sm  a(y-a)f(y)   ^^  _  sin  a(x-a)      T    e^°'(y-^)+  (y)   ^^^ 
J  y'''  a  v/x  ^  y 


a  /x  ^  y  '  a  /x 

'a  -.v 

We  obsei^e  in  the   first  integral  x  ^  y  and  in  the  second  y  S  x.      Hence 
if  we  modify  the  integral  along  the  imaginary  axis  by  replacing  sections  by 
the   above  mentioned  sequence  of  curves   on  which  Im  a  ^  0  we   see  that  the 
singularity  makes  no  contribution. 

III.      Should  a(v)   fail  to  be  even  on  all  or  part  of  the  Im  axis  then  it  is 
necessary  to  add  to  the  discrete  siom  an  integral 

oo 


/         J(x,T)TdT. 


-CD 
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IV.   The  condition  that 

Im  g(v)  <  0 

near  the  negative  imaginary  axis  is  not  necessary.  However^  should  it  be 
positive  then  additional  zeros  of  H  (aa)  will  appear  in  that  region.   Then 
the  proof  would  be  modified  by  considering  a  curve  r„  exactly  analagous  to 
r_  but  near  the  negative  imaginary  axis. 

If  we  replace  the  <  sign  by  a  S  sign  then  zeros  of  H  (aa)  may  occur 
on  the  imaginary  axis.   Near  such  a  point  we  would  modify  the  integral  along 
the  imaginary  axis  by  considering  a  small  semicircle  in  the  right  half  plane 
centered  at  that  point.   This  will  yield  half  the  'ordinary'  residue. 
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Appendix  II 

In  this   appendix  we  will  show  that  the  Watson  seirLes   diverges   for  a 
rather  xinexpected  case.      The  function 

where    p  and  (T  are  different  positive  numbers,    yields   coefficients 


00 


c(v    ) 
^   n 


r  f  (x)H  (x)  ^ 

/         o^         V    ^         X 
^  n 

a 

oo 


such  that  the  seirLes 


oo 


Ec(v   )H     (x) 
.      ^   n     V 
n=l  n 


diverges. 


We  note  that  f  (x)  is  an  analytic  function,  and  that  f  (a)  =  0.  Also, 


for  large  x. 


Hp(x 


1/2 


)-(!.)  -K^-f-j)]  (-^y)  • 


and 


dH  (x) 

P^    ' 

djx 


1/2 


-(^)^  4(^-f-j)](--©) 
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Hence 


dH  \  X/2 


=     ©■ 


JJ^  ' 


Therefore 


11m 
X  — >  00 


df 

-^  -   1   f  (x) 
dx  o 


X  =   0. 


In  spite  of  this,   we  will  show  that 


Ec(v    )H     (x) 
n      V 
n 


diverges . 


We  can  compute  the  coefficients  c(v  )  explicitly  for  f  (x) 


/ 


00  H^(x)H  (x) 
p     V 
n 


dx 


X  H„(x)^  H  (x)  -  H  (x):^  H  (x) 
P^   dx  V        V  ^   dx  p^ 
n       n 


2    2 


X  — >  00 


X  =  a 


■a  H^(a)^  H  (a) 
p^   dx  V  ^ 

n 


2    2 


using  the  above  asymptotic  expressions  and  the  fact  that  H  (a)  =  0.  Also, 

n 


58 


00 


/ 


L  n 


dy 
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V  — >  V 
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n       n 
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2    2 

V   -  V 
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-■  x 
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a 

V 
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00 

/ 

a 

B  (x)B     (x) 

°           d:- 

x 

00 

f 

L    n       -1 

2 
-  dx 

2v 


V^) 


v^  -  p^    r|-B(a)] 


Therefore 


/ 


oo  f  (x)B  (x) 

O^     V 

n 


dx 


2v 


=(vj  = 


n 


^2    2 
(5  -  p 


QD 


/ 


L  n   J  , 
•  dx 


^   V 


H» 


v=v 


,2    2.,  2    2> 
(v^  -  P  )(v^  -  <3-  ) 


In  order  to  estimate  the  c(v  )  and  B  (x),  we  need  the  followi 

n 
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Lemma: 


The   zeros   v     of  H   (a)    lie   on  a  ray 


arg  v^  ~  ^    I  1  - 


1    +    l06 


ae 


Proof: 


For  large   v 


1   sin  jtv  H   (a) 

V 


r(i-v)  r(i+v) 


Equating  and  applying  Stirling's   formula,   ve  find  that 


2v     2v 

-TiV    1 

n       n 

n 

e 

-ae> 

'  2'  'n  2  sin  rt   v 


n 


and 


-2rcos^ 
,a.e.  2rcoB^  -2rTsinY 

(— )        r  e 


-jtvi 


2  sin  nv. 


in  vhlch  we  set 


If 
V  =  re 
n 


This  can  be  satisfied  vhen  cos  f  -is  small.  Furthermore,  we  note  ^  must  be 
close  to  +  rt/2.   We  set 


^■1 


H 


and  find  that 


.  ko  - 


exp  <  2r 


6   log  I  -    (|  -  5) 


or 


6  '- 


2 


2r 

1  +  log  — 

^  ae 


or 


arg  v^  ~  2    [^ 


1    +    l0£ 


2r 
ae  ■ 


We  vlll  now  estimate  H 


V    (x) 


H     (x)   =  ,   .  / 
V   ^  i  sin  Tti 


-jtv  i 
J        (x)        J     (x)e       "" 


-V 


r(i-v  )         r(i+v  ) 

^       n  ^       n 


i   sin  rtv 


-V  -  V 

n  n 

r(i-v^) 


V  V  -JtV    1 

n         n  n 

(f)      (f)        ' 


r(i+v  ) 

^        n 


n     -rtv  1    r 


i  sin  jtv 


n 


/a«.  n 

r(i+v  ) 

^        n 


-V  V 

n  n 

(f)    '(f) 


(since 

H^    (a)   =  0). 
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V 

n 

r(i+v  ) 

n 

n     -V   -1/2 
(3         ^^  ^ 


1  ,8, 

•2'        'n 


1         -1/2  /ae   X    ° 


const  V 


-1/2 


n 


- 1+1 


Tlierefore: 


T.  f        \  ^  -1/2      /Xn       " 

H     (x)  ,->  const   V  (— ) 

n 


We  nov  obtain  estimateB  for 


"Sv      V 


H,(a) 


V  =  V 
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H  (a) 


V  Jtv 
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and 


hE   (a) 

V 


2 


v=v 


^2v  '^ 
n 


^2v  ' 


1  -  log(|^) 


*  -1/2  -, 
const  V     log  V  . 
n        n 


Therefore 


2v  H  (x) 
n   V  ^ 
n 


c(v  )H  (x)  = 
n  V 
n 


Sv"  V 


H  (a) 


2   2 
cr   -p 


/-  2  2,,  2  2, 
(v  -p  (v  -o-  ) 
^  n  '^  ^  n 


v=v 


n 


•^  const 


V  log  V 
n      n 


n 


vhlch  Is  unbounded  for  all  x  >  a.  We  conclude  that  the  series 


^     ^n  \   ^^) 
n 


must  diverge  for  all  x  >  a. 
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